Search for x Cj (2F) from Higher Charmonim El Transitions and X,Y,Z States 
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We calculate the El transition widths of higher vector charmonium states into the spin-triplet 
2P states in three typical potential models, and discuss the possibility to detect these 2P states 
via these El transitions. We attempt to clarify the nature of some recently observed X,Y,Z states 
by comparing them with these 2P charmonium states in these El transitions. In particular, the 
calculated branching ratios of -0(4040), -0(4160) — ¥ Xc,1 (J=0,l,2) are found to be in the range of 
l(T 4 -lCr 3 , and sensitive to the 3S-2D mixing of ^(4040) and ^(4160). The mixing angle may be 
constrained by measuring 0(4040), 0(4160) — > Z(3930)7, if Z(3930) is identified with the Xc 2 state, 
and then be used in measuring Xc ± states. These processes can be studied experimentally at e + e~ 
colliders such as BEPCII/BESIII and CESR/CLEO. 

PACS numbers: 12.39. Jh, 13.20.Gd, 14.40.Pq 
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INTRODUCTION 

Since the discovery of J/ip in 1974 jl], |^], a lot of 
charmonium states had been found in the last cen- 
tury. Among them, the vector states "0(4040), -0(4415), 
and 0(4160), which are commonly assigned as 0(35*), 
ip(4:S) and i{,>(2 3 Di) respectively, can be directly pro- 
duced through e + e~ annihilation into one photon, and 
thus can be readily detected at e + e~ colliders like 
BEPCII/BESIII [g and CESR/CLEO. 

Aside from these vector resonances themselves, it is 
also interesting to detect the products via El transitions 
of these vector resonances into lower charmonium states. 
Especially, the decay channels -0(4040,4160,4415) 
Xcjli J = 0,1,2 can be used to detect the 2P char- 
monia x Cj an d to study the properties of these missing 
states. The importance of experimental establishment of 
these 2P charmonia consists in at least two aspects. On 
one side, the properties of x CJ are important to clarify 
the calculations in various potential models and coupled- 
channel models (see, e.g. || and references therein). On 
the other side, the x c , could be the candidates of some 
of the recently observed charmonium-like states, the so- 
called "X,Y,Z" states (for a review see e.g. M). 

According to potential model estimates, the spin- 
triplet 2P states lie between 3.9 and 4.0 GeV in the 
charmonium family [^|-^). Experimentally, five charmo- 
nium(likc) resonances around 3940 MeV have been found 
recently. The Z(3930) § observed by the Belle Col- 
laboration in 2006 in the 77 fusion experiment with 
a mass 3929 ± 5 ± 2 MeV is identified with the x'c 2 - 
The X(3915), which was also produced in the 77 fu- 
sion experiment Jl0[ and detected in the J/0w channel 
with a mass about 3915 MeV, is possibly the x Co FA - 
The Y(3940) and Y(3915), which were detected in the 
B — > J j ipuiK process by the Belle Collaboration |l2| and 



the BaBar Collaboration (pj separately are another can- 
didates for x Cl - The X(3940), which was found by the 
Belle Collaboration H| in the recoiling spectrum of J/ip 
in the e + e _ annihilation process e + e~ — > J/0 + X and 
e+er -> J/i' + DD*, seems not to be a 2P state fy| H). 
Another well known state, the X(3872), which was first 
found in the J/ip tt + tt~ invariant mass distribution in the 
B meson decay around 3872 MeV with J PC = 1++, 
might be interpreted as the D*°D° + c.c. molecule due to 
the closeness of its mass to the D*°D° threshold. But its 
large production rate in p — p collisions at the Tevatron 
and some properties seem to support that it could be a 
compact bound state, such as the 2P charmonium x Ci , 
or a mixture of the x Cl with the D*°D° + c.c. molecule, 
despite of its lower mass. In fact, the mass of Xci(2P) 
can be lowered to below 3.9 GeV if the color screening 
effects and coupled channel effects are considered |l|. Il5| . 
So it is interesting to examine in the El transitions of 
higher charmonia if the X(3872) can be seen by having 
the x component in its wave function. 

Because of the phenomenological importance of the 
x' c j states mentioned above, we will study the produc- 
tion of these states in the El transitions of higher vec- 
tor charmonia, say, 0(4040), 0(4160) and 0(4415). The 
El transition width can be estimated by potential mod- 
els. Various potential models predict various transition 
widths. However, the transition widths of -0(4415) — s- 
X c , 7 are usually small because of the smallness of the 
overlap integral between the wave functions of 4S and 
2P states. On the other hand, the transition widths 
of 0(4040,4160) — > Xcjli are relatively large (tens to 
hundreds of KeV as those shown in Table [j] and Ta- 
ble O) and the corresponding branching ratios are about 
10 — 10~ 3 . So the x Cj ma Y be detected at the upgraded 
BEPCII/BESIII through these channels. Note that the 
El transition width depends on the phase space which is 
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determined by the mass of x CJ ■ Since the x C2 nas been 
identified with Z(3930), we can compare various poten- 
tial model predictions with the experimental data of x C2 
to see if it can be detected in the El transitions of higher 
charmonium states. This comparison may provide some 
hints in searching for the other two \ c states. 

We introduce three typical potential models in section 
II, and calculate in section III the El transition widths 
of -0(4040,4160,4415) -> x CJ 7 with both lowest- and 
first-order wave functions in the non-relativistic expan- 
sion within these models. We discuss the possibility for 
producing x from El transitions of higher excited char- 
monium states and compare them with those "X,Y,Z" 
states in section IV, where the effects of S-D mixing of 
-0(4040) and -0(4160) on the El transition widths are also 
considered. A summary is given in section V. 



THE POTENTIAL MODELS 

There are many phenomcnologically successful poten- 
tial models in the literature. Among them the Cornell 
model ||(here we mark it by Model II) is well known, 
which describes the charmonium system quite well. How- 
ever, the predicted masses of higher charmonium states 
seem to be larger than their experimental values [(TH). A 
distinct example is the mass of x C2 which is about 40 
MeV larger than the experimental value. The screened 
potential model (see Ref Q and references therein) (here 
we mark it by Model I) was proposed to lower the masses 
of higher charmonia. So we take it here to estimate the 
El transition production of x n , ■ The third model we take 
was proposed by Ding et al. |IJ|, [l?J] (here we mark it by 
Model III), in which the Coulomb potential has a run- 
ning coupling constant. The Hamiltonian in these models 
have the form of H = — ^- + Vy (r) + V~s (r ) , where Vy (r) 
is vector potential and Vs(r) is scalar potential and m c 
is the mass of charm quark. 

The potentials in Model I pl| are: 



Vy(r) 



o » V s( r ) 

3 r 



1 - e _/ir 



(1) 



where fi is the screening factor which makes the long 
range scalar part Vs(r) become flat when r 3> — and 
still linearly rising when r <C — , A is the linear potential 
slope, and ac is the coefficient of the Coulomb potential. 
The model parameters are chosen following Ref Ea|: 



a c = 0.5007, 
fx = 0.0979GeF, 



= 0.21Ge^ 2 , 
1.4045GeV, 



(2) 



where ac ~ a s (m c v c ) is essentially the strong coupling 
constant at the scale m c v c . Here fi is the characteristic 
scale for color screening, and 1/ fi is about 2 /m, implying 
that at distances larger than the static color source 



in the cc system gradually becomes neutralized by the 
produced light quark pair, and string breaking emerges. 
The potentials in Model II are: 

Vv{r) = -^f, V s (r) = Xr, (3) 

with model parameters taken similar to Ref. [Q: 
a c = 0.5461, A = Q.lA2bGeV 2 , m c = 1.4794GeV 
The potentials in Model III ||, arc 
in 1 1 - Ar 



Vy(r) 



V s {r) =Xr-C, (4) 



25 ln(Ar) 1 + Ar ' 
with parameters 

A = 0A7GeV, A = 0.22GeV 2 , 
m c = 1.84GeV, C = -0.975GeK (5) 

The potentials above are used to calculate the lowest- 
order and the first-order non-relativistic wave functions. 
For the first-order relativistic corrections to the wave 
functions, we include the spin-dependent part of H$s, 
Hls, Ht and the spin- independent part Hsi as pertur- 
bations. 

The spin-spin contact hypcrfinc interaction is 



L 7* 7* 2 

HsS = T, — • 52 V Vy. 

Ami 



(G) 



For Coulombic vector potential, V 2 (-) cx <5 3 (r) which 
gives too large a mass splitting of J/tp — r/ c , so we make 
a substitution as in Ref. |7|1 for Model I and Model II: 



Hss = 3 ^ 7rQ |' 7 Mr) S c ■ S E 
9mi 



(7) 



where 6 a (r) = {a/^fe- y2r2 and a = 1.362Gey in 
Model I and a = 1.0946Gey in Model II. 
The spin-orbit term is 



H 



1 



LS 



2m 2 c r 

and the tensor force term is 
1 , 1 „' 



(3V v (r)-V s (r))L-S, (8) 



Ht — 



(-Vy(r)-V v (r))S 12 , 



(9) 



where 5i2 = 3(<t1 • r)(<72 • ?) — d[ ■ d^.- 

The spin-independent part is a bit complicated. We 
take the form as Ref. E0(|: 



Hsi = - 



P 4 



1 



V 2 V v (r) 



4m 3 Am 2 
U{A-W(r)9-P 2 }} 



2m 



2ml 



(10) 
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where Pi and P2 arc momenta of c and c quarks in 
the rest frame of charmonium, respectively, which sat- 
isfy Pi = —P 2 = P, 3 is the unit second-order tensor, 
and {{ }} is the Gromes's notation 

{{A-&-B}} = ^{AB : M+A-VIB+B-VIA+Vl : AB), (11) 

where 3? is a second-order tensor. 

Note that we do not include the contributions from 
the scalar potential in Hsi since it is still unclear how to 
deal with the spin-independent corrections arising from 
the scalar potential theoretically. 

El TRANSITION WIDTHS 

El transitions of higher excited S- and D-wave char- 
monium states are of interest here because they can be 
used to produce and identify P-wave states. For the El 
transition width for charmonium, we use the formula of 
Rcf. |l|: 

r^n^Lj-^'+^+T) 

= ^C fi 5 ss ,eta\{f\r\i)\ 2 E* (12) 

where E 7 is the emitted photon energy. 
The spatial matrix element 

poo 

<f\r\i>= / R f (r)Ri(rydr, (13) 

JQ 

involves the initial and final state radial wave functions, 
and the angular matrix element C/% is 

C fi = max(L, L') (2J' + 1) { j' ^ 1 } ■ ( 14 ) 

We are only interested in initial states with J PC = 
1 , i.e., 0(4040), -0(4415) and *0(416O), since they can 
be easily produced in e + e~ annihilation and can transit 
into spin-triplet 2P states by emitting a photon. 

For the masses of initial and final states used to calcu- 
late E~ in above three models, we take the central values 
from PDG(2010) || if the states are well established 
experimentally The mass of x Co 1S supposed to be 3915 
MeV, while for the mass of x Cl we choose two different 
values: 3872 MeV and 3915 MeV. 

The calculated results with lowest-order wave func- 
tions are listed in Table |. The results of Barnes, et al. M , 
which arc similar to the Model II and the results of God- 
frey, et al. H , with a relativized Cornell model, are also 
listed in Table | for comparison. From Table [| one can 
see that the widths r(-0(4415) — > x' c jl)> are vei T small 
due to large cancelation in the overlap integral between 
wave functions of 4S and 2P states, and therefore are 
very sensitive to the model details. On the other hand, 



the predictions for the widths r(^(4040, 4160) x' c jl) 
are large and insensitive to model details, resulting in 
quite steady values in different models. 

The results obtained with first-order wave functions 
are listed in Table || From (7-10), one can see that the 
corrections to the non-relativistic potential involve some 
derivative terms, which make the potential to be more 
attractive towards the origin. As a result, the wave func- 
tions with relativistic corrections will be thinner than 
those without relativistic corrections. This effect usually 
reduces the spatial matrix elements < f\r\i > defined 
in ([l3]), which can be seen directly by comparison be- 
tween the results of r(0>(4O40, 4160) -> x' c jl) listed in 
Table | and in Table ||. However, the relativistic correc- 
tions can also change the node structures of the wave 
functions of higher exited states, such as -0(4415), and 
make the cancelation in the overlap integral between the 
wave functions of 4S and 2P states more modest, and this 
can be seen in Table [n| On the other hand, the transi- 
tion width is proportional to the factor Ez, which favors 
initial states with higher masses. Thus, the decay widths 
r(^(4415) — > x'cjl) listed in Table |l| become larger. But 
we should mention that these results are not very reliable 
and are more sensitive to the model details than those of 

r(V;(4040,4160)^Xcj7)- 

We calculate the El transition branching ratios of 
0(4040), 0(4160), and 0(4415) with their total width 
taken from PDG(2010) Since the errors of the total 
widths are relatively small for these states, we only take 
the central values of the total widths in calculating the 
branching ratios and do not consider the errors. 

DISCUSSIONS ON XYZ STATES 

In this section, we fucus on the implication of the re- 
sults of Br(0;(4O40, 4160, 4415) -> x' c jl) 011 searching for 
XYZ states in these channels. The 3S-2D mixing effects 
of -0(4040) and -0(4160) are also considered in details. 

Z(3930) 

The Z(3930) was found by the Belle Collaboration g 
in the process 77 — > DD with 

A/(Z(3930)) = 3929 ± 5 ± 2 MeV, (15) 
T(Z(3930)) = 29 ±10 ±2 MeV, (16) 
r 7T 8(Z(3930) -> DD) = 0.18 ± 0.05 ± 0.03 KcV, (17) 

and confirmed by the BaBar Collaboration |2j| with 

M(Z(3930)) = 3926.7 ± 2.7 ± 1.1 MeV, (18) 
r(Z(3930)) = 21.3 ±6.8 ±3.6 MeV, (19) 
r 7T 8(Z(3930) -> DD) = 0.24 ± 0.05 ± 0.04 KcV. (20) 
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TABLE I: El transition widths and branching ratios of charmonium states with the lowest-order wave functions in various 
potential models. The masses and total widths of the initial states used in the calculation are the PDG ]22| central values, while 
the masses of the final states in the Model I-III calculations are denoted by the numbers in the parentheses. For comparison, 
the results of Refs. M and M are also listed. 



Process 


< f\r\ 
— «^ 1 ' 1 


i > {GeV' 1 ) 




k (MeV) 






thy 


(keV) 






Br thy (xlO- 4 ) 




Tntial 


Final 


I 




III 


Ours 


Ref.[[7 


1 Rpf flSfl 


I 


II 


III 


Ref.f? 


Ref.l 

-LL_ 


I 


II 


III 


Rpf 17 

licl. II 1 


1 Ref.il 






-4.9 


A A 


A 


109 


R7 




74 


59 


36 


14 


48 


9.3 


7 A 


A ^ 


1 & 


6.0 




y' r (3872) 


-4.9 


-4.4 


-3.4 


164 


113 


145 


151 


122 


74 


39 


43 


18.9 


15.3 


9.3 


4.9 


5.4 




x' Cl (39i5) 


-4.9 


-4.4 


-3.4 


122 


113 


145 


63 


51 


31 


39 


43 


7.9 


6.4 


3.9 


4.9 


5.4 




xl ( 3915 ) 


-4.9 


-4.4 


-3.4 


122 


184 


180 


21 


17 


10 


54 


22 


2.6 


2.1 


1.3 


6.8 


2.8 


0(4160) 


Xc 2 (3929) 


5.0 


4.6 


3.8 


218 


183 


210 


12 


10 


7.3 


5.9 


6.3 


1.2 


0.97 


0.71 


0.57 


0.61 




Xc!(3872) 


5.0 


4.6 


3.8 


271 


227 


234 


355 


299 


210 


168 


114 


34.5 


29.0 


20.4 


16.3 


11.1 




Xl 1 (3915) 


5.0 


4.6 


3.8 


231 


227 


234 


219 


185 


132 


168 


114 


21.3 


18.0 


12.8 


16.3 


11.1 




xl (3915) 


5.0 


4.6 


3.8 


231 


296 


269 


292 


247 


173 


483 


191 


28.3 


24.0 


16.8 


46.9 


18.5 


0(4415) 


Xc 2 (3929) 


-0.013 


0.093 


-0.028 


465 


421 


446 


0.04 


2.1 


0.2 


0.62 


15 


0.006 


0.34 


0.03 


0.1 


2.4 




x' Cl (3872) 


-0.013 


0.093 


-0.028 


515 


423 


469 


0.04 


1.7 


0.2 


0.49 


0.92 


0.006 


0.27 


0.03 


0.08 


0.15 




Xc!(3915) 


-0.013 


0.093 


-0.028 


477 


423 


469 


0.03 


1.4 


0.1 


0.49 


0.92 


0.005 


0.23 


0.02 


0.08 


0.15 




xl (3915) 


-0.013 


0.093 


-0.028 


477 


527 


502 


0.01 


0.45 


0.04 


0.24 


0.39 


0.002 


0.07 


0.006 


0.04 


0.06 



TABLE II: El transition widths and branching ratios of charmonium states in various potential models with the first-order 
wave functions. The masses and total widths of the initial states used in the calculation are the PDG j2^] central values, while 
the masses of the final states in the Model I-III calculations are denoted by the numbers in the parentheses. 



Process 


</ 


r\i > (GeV- 1 ) 


k(MeV) 


J 


r thy (keV) 




B 


r thy (xlO- 4 ) 


Initial 


Final 


I 


II 


III 




II 


III' 


I 


II 


III 


0(4040) 


Xc 2 (3929) 


-4.3 


-3.9 


-3.1 


109 


56 


47 


29 


7.0 


5.9 


3.6 




xk(3872) 


-3.7 


-3.4 


-2.8 


164 


88 


72 


50 


11.0 


9.0 


6.3 




Xc,(3915) 


-3.7 


-3.4 


-2.8 


122 


37 


30 


21 


4.6 


3.8 


2.6 




Xc n (3915) 


-3.0 


-2.7 


-2.5 


122 


7.9 


6.2 


5.3 


0.99 


0.78 


0.66 


0(4160) 


Xc 2 (3929) 


4.3 


4.1 


3.4 


218 


9.2 


8.2 


5.9 


0.89 


0.80 


0.57 




x' ci (3872) 


3.6 


3.4 


3.2 


271 


189 


169 


147 


18.3 


16.4 


14.3 




x' cl (3915) 


3.6 


3.4 


3.2 


231 


117 


105 


91 


11.4 


10.2 


8.8 




Xc (3915) 


2.7 


2.6 


2.9 


231 


89 


81 


97 


8.6 


7.9 


9.4 


0(4415) 


Xc 2 (3929) 


-0.42 


-0.13 


-0.19 


465 


42 


4.1 


8.7 


6.8 


0.66 


1.4 




x' ci (3872) 


-1.1 


-0.77 


-0.39 


515 


219 


116 


30 


35.3 


18.7 


4.8 




Xci(3915) 


-1.1 


-0.77 


-0.39 


477 


174 


92 


24 


28.1 


14.8 


3.9 




Xc (3915) 


-1.8 


-1.4 


-0.65 


477 


164 


105 


22 


26.5 


16.9 


3.5 



The production rate and the angular distribution in 
the 77 center-of-mass frame suggest that this state is the 
previously unobserved x C2 || Its mass, however, is 
about 40-50 MeV larger than the commonly predicted 
value in the quenched potential model (see, e.g. M). 
A lower mass can be obtained by considering the color 
screening effect described in Model I jL5| in which the 
predicted mass is 3937 MeV. 

Since Z(3930) is established as the candidate of 
X C2 , searching for Z(3930) in the El transitions of 
■0(4040, 4160, 4415) is important to further verify this as- 
signment and can also be a good criterion in searching 
for and identifying other x' c j states in these transitions. 



From Table § and Table [fj, we can see the transition 
width of 0(4040) -> x' C2 7 is 36 " 74 KeV wit h the lowest- 
order wave functions and 29-56 KeV with the first-order 
wave functions in our calculations within Models I-III 
and in Ref. The corresponding branching ratio is 4.5- 
9.3 xlO -4 and 3.6-7.0 x 10 -4 , respectively. The branching 
ratio of order of 10~ 4 is encouraging to detect x C2 m 
■0(4040) El transitions. Note that the results of Ref. @] 
is notably small than our results. This is mainly because, 
the mass of x C2 use d in Ref. @ is larger than ours and 
the corresponding energy of the emitted photon is much 
smaller than ours. 

The branching ratio of ^(4160) — > x C2 7 i s about one- 
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fifth of that of ^(4040) x C2 7 but is still close to 
1 x 1CT 4 . The branching ratio of -0(4415) — > x C2 7 1S 
sensitive to the model details. So it is difficulty to pre- 
dict how large is the branching ratio of ■0(4415) — > x C2 7 
in potential models. 



X(3915),Y(3940),Y(3915) 

The X(3915) (lO), Y(3940) (TJ and Y(3915) [|TJ arc 
all observed in the invariance mass distribution of J/-0 w 
in processes 

77^ X(3915) -^J/0w, (21) 

B -> i^F(3940), Y(3940) -+ J/ipw, (22) 

B ^ KY (3915), Y(3915) J/ipu, (23) 



with masses 



M(X(3915)) : 
M(y(3940)) : 
M(y(3915)) : 

total widths 

r(X(3915)) 
T(Y(3940)) 
T(Y(3915)) 

and partial widths 



3915 ± 3 ± 2 MeV, 
3943 ±11 ±13 MeV, 
3914.6±|;1 ± 2.0 MeV, 



17 ± 10 ± 3 MeV, 
87 ± 22 ± 26 MeV, 
34±* 2 ± 5 MeV, 



(24) 
(25) 
(26) 



(27) 
(28) 
(29) 



T 77 x B(X(3915) - J/* U ) ={l\± ^ J f P " £ helidty . 2 (30) 

23(5 ->• AT y (3940)) x 23(Y(3940) -> J/^w) = (7.1 ± 1.3 ± 3.1) x 10" 5 , (31) 

B{B+ K + y(3915)) x B(Y (3915) -> J/0 w) = (3.5 ± 0.2 ± 0.4) x 10~ 4 , (32) 

6(5° if y(3915)) x S(y(3915) -> J/0w) = (3.1 ± 0.6 ± 0.3) x 10~ 4 . (33) 



Although the differences of masses and total widths 
of the three signals are within 2a, especially those of 
X(3915) and Y(3915) are less than one a, it is not clear 
whether these three signals come from the same particle. 
BaBar @ considers Y(3915) and Y(3940) as the same 
state since the smaller values of both the mass and to- 
tal width of Y(3915) derived from fitting data by BaBar 
can partially be attributed to larger data sample used 
by BaBar, which enable them to use smaller J /ip to mass 
bin in their analysis [p4f . Y(3940) and X(3915) are also 
considered to be the same state in Rcfs. |1C], |24|-|2q1. 

Besides whether they are the same particle or not, 
there are no decisive interpretations of these states. That 
they are detected in the J/0 uj channel but not in the 
DD or DD* channel makes people suspect they are 
not conventional charmonium states. Ref. jlO| argues 
that X(3915) is not an excited charmonium state but 
favors the prediction of D*D* bound state model p7) . 
The Y(3915) and Y(3940) arc also interpreted as D*D* 
molecular states by |28|, |29| . 

However, Liu ct al. |1J] pointed out that the node effect 
of excited wave functions may change the open charm de- 
cay widths dramatically. They calculated the open charm 
decay of x c , an d Xc, m the 3 Pq model, and argued that 
X(3915) is the x co state. 

If they are conventional charmonia, they are probably 
the x Cj states since their masses are consistent with po- 



tential model predictions and they have positive 

charge parity since they decay to J/ipoj. 

The X(3915), which is observed in 77 fusion, may be 
X ca or x C2 ■ -But the mass of X(3915) is about 2a different 
from that of Z(3930) and from (JTj) and ©. Further- 
more, if X(3915)is the same meson as Z(3930), we would 
get B(x' C2 -> J/4>uj)/B{x C2 ) -> DD) « 0.1 ±0.04, which 
seems to be too large for charmonium. So X(3915) is un- 
likely to be Xc 2 an d we tend to regard it as a candidate 

fOT x'c - 

The Y(3915), which is produced in B decays, has so 
close mass and total width to the X(3915) that we sus- 
pect they are the same state. However, if Y(3915) is 
not X(3915), then it may be x Cl - Since the El transition 
rate from 0(4040) to x Cl is three to four times larger than 
that to Xcn w ith both lowest-order and first-order wave 
functions if they have the same mass, we can distinguish 
between them by measuring these El transitions. 

The Y(3940) is unlikely to be x' Cl - If it is X Cl , its 
main decay mode ought to be DD* . However, Belle gives 
23(y(3940) -> J/0 w)/B(Y (3940) -> D°D*°) > 0.71 at 
90% CL p0[ , which disfavor the x' ci assignment. But 
Y(3940) might still be a candidate for % . 

We may detect and identify them in the El transi- 
tions of higher charmonium states. We can see from 
Tabic H and Table || that the transition widths are 
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0(10) KeV for 0(4040) -> Xco i 7 and 0(100) KcV for 
■0(4160) — > Xc 1 7' corresponding to branching ratios of 
order 10~ 4 and 10~ 3 , in all listed potential models. The 
effect of corrections of wave functions is moderate for 
-0(4040) and -0(4160) transitions. The large branching 
ratios may make the detection of x Co i possible. 

The transitions of -0(4415) — > \c 1 7 are model de- 
pendent and sensitive to the node structure, just like 
0(4415) — > x C2 7 as we nave remarked. 

X(3872) 

The X(3872) was first observed by Belle in 
the J/0 7r + 7r~ invariant mass distribution in B + — > 
K + J/t/j tt + it~ decay as a very narrow peak (Fx < 2.3 
MeV) around 3872 MeV. The mass of A (3872) in the 
J/0 ir + ir~ mode was recently updated by CDF Collab- 
oration ||l| as 

M(X(3872)) = 3871.61 ± 0.16 ± 0.19 MeV, (34) 

which is very close to the D°D*° threshold m(D°D* ) = 
3871.81 ± 0.36 MeV Q. Moreover, analyzes both by 
Belle |3j| and CDF ]34| favor the quantum number 
J PC = 1 ++ . The mass seems to be too small for a 
J PC = 1 ++ Xci charmonium, but the color-screening ef- 
fects and the coupled channel effects may lower its mass 
towards 3872 Mev |l5| , as has been declared in Section 
I. 

There are lots of possible explanations for 
X(3872)(sec jij, H) for a review). Aside from the 

the D°D*° molecular state, 
or a mixed 1 ++ 




most popular one, i.e. 
the 1 ++ charmonium 

charmonium- Z)°Z)* state [|37[ |3§|] for X(3872) was also 
proposed. More data samples are needed to distinguish 
between various explanations. 

Recently, an analysis of the uj — > 7r + 7r~7r° spectrum in 
the decay B KX KJ/tpuj performed by BaBar |||] 
claimed that the J PC of A(3872) might disfavor 1 ++ , as 
had widely been accepted, but favor 2 If this result 
is confirmed, the natural assignment is the 1 J>2 char- 
monium i] C 2- However, the mass of this D-wavc state is 
about 3.80-3.84 GeV in the potential models which 
is too low to be the candidate of A (3872). Besides, some 
recent theoretical studies on the properties of r] C 2 indicate 
that it is not apt to the the candidate of A(3872) [|o|-|3). 
Therefore, we will ignore this possibility in the following 
analysis. 



We are interested here in detecting X(3872) in the El 
transitions of higher charmonia, if X(3872) is the 2P char- 
monium Xci or contains some 2P charmonium component 
in its wave function. One can see from Table || that the 
branching ratio for 0(4040) x' Cl 7 is (0.93-1.89) x 10" 3 

and that for 0(4160) -> Xci 7 is ( 2 - 04 ~ 3 - 05 ) x 10_3 in 
our calculation with zero-order wave functions. The large 
branching ratios may enable us to find x Cl m the e + e~ 
machines and compare it with X(3872). Note that Rcf. Q| 
and Ref. Q give smaller branching ratios. It is partly be- 
cause they take a larger mass for x in calculations. If 
they take the mass of x Cl to be the same as us, the dif- 
ferences will diminish. This means the branching ratios 
are not sensitive to models. The calculated results with 
relativistically corrected wave functions are a bit smaller 
but still quite large (see Table ||). It indicates that the 
results are not sensitive to the nodes of wave functions 
and should be reliable. 



3S-2D Mixing 

The 0(4040) and 0(4160) are commonly assigned as 
0(3 3 Si) and 0(2 3 £>i) respectively. Therefore, the above 
results of El transitions are all based on these sim- 
ple assignments. However, the observed leptonic width 
T ee (4040) w r ee (4160) is inconsistent with this picture. 
The simplest explanation is that they are roughly 1:1 
mixtures of 0(3 3 Si) and 0(2 3 Z?i). Neither the ten- 
sor force nor the coupled channel effects can cause such 
strong mixing (see Ref. [0 and references therein) and 
the mixing mechanism remains unknown. Here we do not 
consider the mixing mechanism, and simply assume that 
they are mixtures of 0(3 3 Si) and 0(2 3 J0) with a mix- 
ing angle 8, and calculate how the El transition widths 
vary with 9. In this case, the 0(4040) and 0(4160) are 
expressed as 



1^(4040)) = 
10(4160)) = 



|3 3 Si)cos6> + 
-|3 3 Si)sin6» 



|2 3 Di)sin0, 
f |2 3 Di)cos6». 



(35) 
(36) 



Using the data of leptonic decay widths of 0(4040) and 
0(4160) p3] as inputs, one can determine the mixing 
angle 9 in the potential models. It is about —35° or 
+55° in all the three models used in Sec. II. And this is 
consistent with earlier estimates of the mixing angle 45 



The corresponding El transition widths parameterized 
by the mixing angle are 
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r(0(4O4O) -> Xco 7) 



— ae 2 c k 3 (cos 2 9 < 2 3 P |r|3 3 ,S'i > 2 -2^/2cos6sin6 < 2 3 P |r|3 3 S"i >< 2 3 P |r|2 3 P>i > 
+2sin 2 9 < 2 3 P |r|2 3 P> 1 > 2 ), 



(37) 



r(0(4O4O) -> Xci7 ) 



r(0(4O4O) Xc2 7) 



-ae 2 c k 3 (cos 2 6 < 2 3 Pi|r|3 3 S*i > 2 +V2cos9sin6 < 2 3 Pi|r|3 3 S'i >< 2 3 Pi|r|2 3 P>i > 
+ ^sin 2 6 < 2 3 P 1 \r\2 3 D 1 > 2 ), 



20 
27 



ae 2 c k 3 {cos 2 9 < 2 3 P 2 |r|3 3 S'i > 2 -—cos6sin6 < 2 3 P 2 |r|3 3 Si >< 2 3 P 2 |r|2 3 P» 1 > 



+ -sm 2 fl < 2 3 P 2 |r|2 3 J D 1 > 2 ), 
50 



(38) 



(39) 



r(0(416O) -> Xco 7) 



r(V>(4160) -> Xci 7) 



r(0(416O) -> Xc2 7) 



ae 2 fc 3 (2cos 2 6> < 2 3 P |r|2 3 P>i > 2 +2V2cos0sin6 < 2 3 P |r|2 3 P>! >< 2 3 P |r|3 3 Si > 



27 

+sin 2 6 < 2 3 P |r|3 3 5i > 2 , 

^-ae 2 c k 3 {]-cos 2 9 < 2 3 P 1 \r\2 3 D 1 > 2 -V2cos6sin6 < 2 3 P 1 \r\2 3 D 1 >< 2 3 Pi|r|3 3 S"i > 
9 2 

+sin 2 6 < 2 3 P 1 |r|3 3 S'i > 2 ), 



(40) 



(41) 



20 1 \p2i 

T^^ 2 c k 3 {—cos 2 e < 2 3 P 2 \r\2 3 D 1 > 2 +^—cos9sin0 < 2 3 P 2 \r\2 3 D 1 >< 2 3 P 2 |r|3 3 5i > 
27 ' 50 5 

+sin 2 e < 2 3 P 2 |r|3 3 5i > 2 ), 



(42) 



r 



We use the lowest-order wave functions calculated in 
model I-III and impose m x > = 3929 MeV, m^j = 
3872 MeV and m< = 3915 MeV. The results are similar 

Acq 

in the three models, i.e., the El transition widths reach 
their maximum or minimum values almost at the same 
mixing angle in the three models. As an example, we 
display the results in Model I and Model II in Figs. 
||, and do not show the similar results in Model III for 
simplicity. 

We see that the decay width of V(4040) -> xL, (3929)7 
reaches its maximum of about 75 KeV corresponding to a 
branching ratio of 9.3 x 10~ 4 near = 10°. While the de- 
cay width of •0(4160) — > x C2 (3929)7 reaches its minimum 
(zero) near 10° and its maximum of about 600 KeV near 
—80° at which •0(4160) is almost a pure 0(3 3 S"i) and the 
corresponding branching ratio is about 6 x 10~ 3 . 

The decay width of V(4040) ->• (3872)7 reaches its 
maximum of about 250 KeV corresponding to a branch- 
ing ratio of 3.1 x 10~ 3 near —35°, which is about 2 
times larger than that of non-mixing, and reaches its 
minimum (zero) near 55°. While the decay width of 
0(4160) — > Xci (3872)7 reaches its minimum (zero) near 
—35° and reaches its maximum of about 1050 KeV cor- 
responding to a branching ratio of about 1% near 55°. It 
is very interesting to note that the two values happen to 
be corresponding to the mixing angles determined by lep- 
tonic decay widths of ^(4040) and ^(4160). That means, 



one of the two channels must be enhanced by the 3S-2D 
mixing if the mixing mechanism just affect the leptonic 
decays and the El transitions by simply mixing the wave 
functions. 

The decay width of 0(4040) -> x' co (3915)7 reaches its 
maximum of about 52 KeV corresponding to a branch- 
ing ratio of 6.5 x 10~ 4 near 55°, which is about 2.5 
times larger than that of non-mixing and reaches its 
minimum (zero) near —35°. While the decay width of 
0(4160) — > x c (3915)7 reaches its minimum (zero) near 
55° and reaches its maximum of about 450 KeV corre- 
sponding to a branching ratio of about 4.5 x 10 ~ 3 near 
-35°. 

Since the decay width of 0(4160) -> \ C2 (3929)7 
reaches its minimum (zero) and 0(4040) — > \ C2 (3929)7 
reaches its maximum near 10°, which is not far from the 
non-mixing case, in our model calculations, we may use 
these two channels to check whether there is substantial 
3S-2D mixing between 0(4040) and ^(4160). Ifxl 2 (3929) 
can be detected in the El transitions from •0(4040) but 
not from 0(4160), then the mixing angle should be small. 
In general, since the Z(3930) is identified with x C2 , the ob- 
served El transition rates to Z(3930) from 0(4040) and 
0(4160) will be useful to constrain the value of 3S-2D 
mixing angle by comparing the measurements and the 
theoretical predictions. When we have a better control 
over the value of mixing angle, we will be in a position 
to study the properties of X(3872) and X(3915), and to 
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FIG. 1: El transition widths of 0(4040) varying with 3S- 
2D mixing angle in Model I. A: 0(4040) -> X l 2 (3929)7, B: 
0(4040) -J- Xcx (3872)7, C: 0(4040) -> xl„(3915) 7 . 



FIG. 3: El transition widths of 0(4160) varying with 3S- 
2D mixing angle in Model I. A: 0(4160) -¥ x'c 2 (3929)7, B: 
0(4160) -J- Xc! (3872)7, C: 0(4160) -> X l„(3915) 7 . 




> 
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Mixing Angle 




-20 20 

Mixing Angle 



FIG. 2: El transition widths of 0(4040) varying with 3S- 
2D mixing angle in Model II. A: 0(4040) -> x'c 2 (3929)7, B: 
0(4040) -> (3872) 7 , C: 0(4040) -S- xl (3915) 7 . 



FIG. 4: El transition widths of 0(4160) varying with 3S- 
2D mixing angle in Model II. A: 0(4160) ->■ xl 2 (3929)7, B: 
0(4160) -S- Xc x (3872)7, C: 0(4160) -> Xc (3915)t 



see if they can be respectively the x Cl (or partially be) 
and Xc ^ comparing the observed El transition rates 
with theoretical predictions for ■0(4040) and -0(4160) de- 
cays to Xci 7- However, we must keep in mind that the 
3S-2D mixing model for 0(4040) and 0(4160) is only a 
simplification for the real situation, since the mixing with 
4S state and the charm meson pairs are all neglected in 
the 3S-2D mixing model. Nevertheless, with this simple 
model we hope some useful information can be obtained 



for these X,Y,Z states by measuring the El transition 
rates of 0(4040) and -0(4160) to x cj (J=2,l,0) charmo- 
nium states. 

We also calculate the El transition widths of 
■0(4040,4160) — > Xcj(l-P)7 in order to see whether these 
can be helpful in determining the 3S-2D mixing angle. 



The results are listed in Table HI. Although the calcu 



latcd widths arc cxpcctcdly small, the transition branch- 
ing ratio of \ Cl 2 — > J /if; 7 are relatively large (about 20% 
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for Xc 2 an d 36% for Xci)- So hopefully it is possible to 
measure them in experiment if the data samples are large 
enough. A special case is '0(4160) — > Xci7 with a branch- 
ing ratio of order 10 ~ 4 , which is quite large, and might 
be easier to be detected. Again, if we have determined 
the mixing angle, it would help us searching for x C; an d 
identify the X,Y,Z resonances. 

SUMMARY 

In this paper, we calculate the El transition widths 
and branching ratios of 0(4040,4160,4415) — > x C; 7 m 
three typical potential models with both lowest- and 
first-order relativistically corrected wave functions. We 
find the transition widths of 0(4040,4160) — > \ CJ are 
model-insensitive and relatively large (tens to hundreds 
of KeV) and the corresponding branching ratios are 
of order 10~ 4 — 10~ 3 , which make the search for \ CJ 
possible at e + e~ colliders such as BEPCII/BESIII and 
CESR/CLEO. This may help us identify some of the re- 
cently discovered X,Y,Z particles, especially the Y(3915), 
Y(3915), Y(3940) and Y(3872). The possible 3S-2D 
mixing of 0(4040) and 0(4160) and its effect on the 
transition widths of 0(4040,4160) — s- x Cl l arG consid- 
ered and found to be important. We find the transitions 
■0(4040,4160) — > Xc 2 7 can ^ c use d to examine whether 
the mixing exists and to further possibly constrain the 
mixing angle; and the transitions 0(4040,4160) — > x Cl l 
and 0(4040,4160) — > X CQ 1 can be used to estimate how 
far is the mixing angle from —35° and 55°, which are de- 
termined in a simple 3S-2D mixing model by the observed 
lcptonic decay widths of 0(4040) and 0(4160). The tran- 
sitions of 0(4040,4160) — > Xcj (1-P)7 are a l so discussed, 
and the 0(4160) — > Xc x l transition is found to have a 
branching ratio of order of 10 -4 , which may be relatively 
easier to be detected. 
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